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The zero temperature core-level photoemission spectrum of a Hubbard system is studied across 
the metal to Mott insulator transition using dynamical mean-field theory and Wilson's numerical 
renormalization group. An asymmetric power-law divergence is obtained in the metallic phase 
with an exponent a{U, Q) — 1 which depends on the strength of both the Hubbard interaction U 
and the core-hole potential Q. For Q < Uc/2, a decreases with increasing U and vanishes at the 
transition {U — > Uc) leading to a symmetric peak in the insulating phase. For Q > Uc/^, a remains 
finite close to the transition, but the integrated intensity of the power-law vanishes and there is no 
associated peak in the insulator. The weight and position of the remaining peaks in the spectra can 
be understood within a molecular orbital approach. 
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I. INTRODUCTION 

When an incident x-ray photon ejects an electron from 
a core-level in a metal, the conduction band electrons feel 
a local attractive potential due to the created hole. It 
was discovered by Anderson^ that the electronic ground 
states before and after the creation of the hole are orthog- 
onal to each other. This many-body effect has dramatic 
consequences in x-ray photoemission spectroscopy (XPS) 
experiments where an asymmetric power-law divergence 
is observedi^i^ For a non-interacting metal, the expo- 
nent of the power-law and the relative intensity of the 
peaks in the XPS spectra are well understood. However, 
the behavior of the power-law divergence in a strongly 
interacting metal has received little theoretical attention 
besides the one-dimensional case<^i^ 

Recently, there have been several XPS studies of 
strongly correlated transition-metal oxides fi^i^ii^ which 
addressed the changes in the core-level spectrum across 
the metal to Mott insulator transition (MIT). It was 
observed in particular that the (strongly renormalized) 
low-energy quasiparticles present in the metallic phase 
close to the MIT strongly modify the satellite peaks 
measured in XPS, an effect also discussed theoretically 
by Kim et al^ In this paper, we provide a detailed 
theory of the core-level photoemission lineshape across 
the metal-insulator transition. A power-law behavior is 
found throughout the metallic phase up to the MIT where 
it is destroyed. Both the exponent and the intensity of 
the power-law are strongly renormalized by interactions 
and two different regimes can be identified depending on 
the ratio between the valence band interaction and the 
core-hole potential intensities. 

The rest of the paper is organized as follows. In Sec. 
|TT]we describe the model that we use in our calculations. 
The numerical results for the XPS spectra across the MIT 
transition are presented and discussed in Sec. IIIII Sec- 
tion llVl contains the results of a molecular orbital approx- 
imation. Finally, we state the conclusions of our study 



in Sec. El 



II. THE MODEL 



The simplest model to study the correlation induced 
MIT is the Hubbard model 
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where cj^ creates an electron with spin a at site i and 
nia = c\^Cicr. A key quantity for this problem is the 
spectral function A{u}) — —-hn[Gii{ijj -t- iO)] that con- 
tains information about the valence-band single particle 
photoemission spectrum. It can be obtained in the limit 
of large lattice coordination using dynamical mean-field 
theory (DMFT)." In the DMFT framework the Hubbard 
model reduces to an Anderson impurity model coupled 
to a non-interacting electron bath that is calculated in a 
self-consistent way. 

In the XPS experiment a photon excites a core elec- 
tron out of the sample and the resulting core-hole inter- 
acts with the band electrons attractively. In the sudden 
approximation considered below, the self-consistent bath 
remains unchanged while the energy of the effective An- 
derson impurity is shifted. The core-hole potential is 
taken as local and momentum independent^i^ 
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where h''{h) creates (destroys) a core-hole at site 0, and 
eti is the core-level energy. In the sudden approxima- 
tion the XPS spectrum is given by the core-hole spectral 
density 
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where the final states with corresponding energies 

{El} satisfy the sum rule J2i,f l(^^/|Oi)P = 1 and [0^) is 
the initial ground state with energy Eq. The spectrum is 
zero for energies larger than the threshold energy ojt = 
Eq — Eq — Ch where Eq is the ground state energy in the 
final configuration. 

In what follows we will consider a Bethe lattice geom- 
etry where the free density of states is given hy^ 




we will set the half-bandwidth D = 2t a.s the energy unit, 
and focus on half-filling {ed — —fJ. — U/2). 



III. NUMERICAL RESULTS 

In the non-interacting case {U 0) the spectral func- 
tion A{uj) is simply is given by Eq. ([3]) and the XPS spec- 
trum then presents a Doniach-Sunjic power-law behavior 
at the threshold^i^iiii^ii^ 

Here Co 1, a = 2(6/7:)^ = 2/7r2 arctan2(2Q), where 
5 = 5{eF) is the phase-shift at the Fermi level due to 
scattering with the core-hole, and the factor 2 is due to 
spin degeneracy. For Q > D/2 a. bound state appears 
in the final configuration that is refiected in the XPS 
spectra as an additional peak at a; — lut ~ —Q- These 
analytical results have been nicely reproduced numeri- 
cally using Wilson's numerical renormalization group 
(NRG) by Lfbero et al}^ Here we will use a different ap- 
proach within the NRG scheme. To improve the accuracy 
at high energies we use a density matrix formulationji^ 
We have 

= ^ S{u;-El+El+e,)J2{i.fW,){i.,\pW'^{i.',Wf), 

where p is the density matrix. To evaluate the matrix 
elements {vf\vi) and {yi\p\v'^ two simultaneous runs of 
NRG are performed (with and without the core-hole en- 
ergy shift). The {vf\vi) are calculated on the first itera- 
tion and then transformed to the new basis at each NRG 
iteration. The spectra are constructed following Ref. [H. 
The high numerical resolution of the NRG at low energies 
allows us to recover the power-law behavior with param- 
eters that reproduce the analytical results (see Fig.[T]). A 
bound state emerges in the final state for Q > D/2 and 
the associated peak in the XPS spectra can be observed 

at W — LOT ^ —Q- 

For a non-zero Hubbard interaction, the spectral func- 
tion A{u!) needs to be calculated in a self-consistent way. 
In what follows we will use the density matrix numeri- 
cal renormalization group (NRG) schem o^^i^^i-*^^ to solve 




FIG. 1: (Color online) XPS spectra in the non-interacting 
case at electron-hole symmetry for different values of the core- 
hole potential Q. Inset: numerically obtained a (dots) and 
analytical expression 2/7r'^ arctan^(2(5) (solid line). 



the DMFT equationsi^iiS and to calculate the XPS spec- 
^ J, 0^13, 20 different values of U. The ratio U/D can be 
varied in some materials across the MIT using chemi- 
cal or mechanical pressure. The ratio U/Q, however, is 
expected to remain approximately constant for a given 
material. In Ref. [1 the spectra in Ru 3d systems were 
fitted using Q = 1.25U. 

Figure shows numerical results for the valence 
band spectral function. A small interaction in the band 
{U <C D) produces only a slight modification of A{uj) 
that close to the Fermi level is approximately given 
by Eq. ([3]). For intermediate values of the interaction 
(1.5 < J7 < J7c — 2.95) the system is in the strongly corre- 
lated regime and the valence band has a three three peak 
structure. The central quasiparticle peak has a width 
~ zD, where z^^ = 1 — 95R[E(cj)]/9w|a;=o and S(a;) is 
the self-energy. The incoherent upper and lower Hubbard 
bands have a width ~ D and their positions are deter- 
mined by the atomic energies. The Fermi liquid behavior 
is restricted to a narrow region close to the Fermi level 
inside the quasiparticle peak which narrows with increas- 
ing interaction as z ~ 0.3{Uc — U). For values of U larger 
than the critical interaction Uc the system is in the Mott 
insulator phase, there is no quasiparticle peak and the 
charge excitations are gapped. 

Within the DMFT framework the interactions in 
the band produce a twofold modification of the non- 
interacting x-ray problem. The atomic level coupled to 
the core-hole potential has now a Hubbard interaction U 
and is hybridized with a renormalized electronic bath. 
To lowest order in U we may neglect the changes in 
the valence band and perform a Hartree-Fock approxi- 
mation to treat the local interaction. We are left with a 
non-interacting x-ray problem with a renormalized core- 
hole potential Q' = Q — {n — l)U/2, where ri > 1 is 
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FIG. 2; (Color online) a) Valence spectra across the Mott 
transition in the electron-hole symmetric case. Note that for 
U = 3.0 the spectrum is gapped while for U = 2.922 there is 
a narrow peak at the Fermi level, b) Core-hole spectra for a 
core-hole potential Q = 1.2 across the MIT showing a power- 
law behavior close to the threshold in the metallic phase, c) 
Same as b) with Q — 6.4. d) and e) Detail of the high energy 
peaks of b) and c) , respectively. 



the charge on the impurity in the presence of the core- 
hole potential. In general we have Q' < Q and we 
therefore expect the bound state peak (if present) to be 
shifted toward the threshold (w — ut ~ —Q') and the 
power-law to become more divergent as a is reduced to 
aHF ~ 2/7r2 arctan2(2Q')- 

In Fig. [2h-e we show XPS spectra for different values 
of the interaction U and the core- hole potential Q. We 
obtain a power-law behavior (straight lines) throughout 
the metallic phase that disappears at the MIT. For small 
U = 0.2 the XPS spectrum can be understood within 
the Hartree-Fock scheme. This approximation, however, 
breaks down rapidly with increasing U as the system ap- 
proaches the strongly correlated regime. We may argue 
in this regime that the width of the quasiparticle peak 
zD acts as an effective bandwidth for the x-ray prob- 
lem. While this allows us to understand the decrease of 
the onset energy for the power-law behavior (see Figs.Oj 
and [2t) it does not give the complete picture as we will 
see below. 

Close to the MIT there are some important differences 



in the spectra for the two values of Q shown in Fig. [2l 
While the power-laws in Fig. Ha for Q = 1.2 < Uc/2 
have a significant increase in the slope when approaching 
the MIT (we have a ^ for f/ ^ C/^), those for Q = 
6.4 in Fig. [2}; have nearly parallel displacements with 
only a very small increase in the slope throughout the 
metallic phase. For Q = 1.2 the high energy peaks shift 
toward the edge with increasing U and disappear at the 
MIT leaving a delta-function peak at the threshold in the 
insulator. For Q = 6.4 > t/c/2 the main difference in the 
XPS spectra just below {U = 2.922) and above {U = 3.0) 
the MIT is the presence of the power-law in the metallic 
phase. The high energy peaks are essentially unchanged 
and the integrated intensity of the power-law vanishes at 
the MIT with a^O. 

Some insight on the origin of these two distinct behav- 
iors can be obtained by analyzing the atomic limit in the 
insulator. For values of Q such that Q < Uc/2, we have 
Q < U/2 throughout the insulating phase, which means 
that the core-hole potential is not strong enough (in the 
electron-hole symmetric case) to overcome the local re- 
pulsion and change the charge locally in the atomic level. 
The initial and final ground states will be identical and 
we therefore expect a peak at the threshold. In the op- 
posite situation, Q > Uc/2, there will be a region close 
to the MIT in which the initial and final states have a 
different occupation of the atomic level. For these values 
of Q, the initial and final ground states are orthogonal 
and no peak is expected at the threshold. 
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FIG. 3: (Color online) Left panel: Exponent a as a function of 
U, for different values of g = 12.8, 2.4, 1.6, 1.2, 0.8, and 0.4 
(from top to bottom). Results from NRG (symbols), Hartree- 
Fock results (solid lines), and parabolic fits f{U) = AQ^{Uc — 
U)^ with A ~ 0.04 (dashed lines). Right panel: integrated 
intensity of the power-law peak for Q = 12.8, 1.6, 1.2, and 
0.4 (from bottom to top). The figure illustrates that the peak 
at threshold disappears at the MIT in two different manners, 
for Q > Uc/2 and Q < Uc/2 (see text). 

The results for the exponent and the integrated in- 
tensity of the power-law are summarized in Fig. [3l The 
exponent a increases with Q and decreases with increas- 
ing interaction. For small U Q,D the main effect of 
the interaction is a compensation (screening) of the core- 
hole potential that can be explained within a Hartree- 
Fock approach (soHd Hues in Fig. [3]). For large Q » Uc, 
a is close to its maximum value 1/2 and the interaction 
produces only a small reduction throughout the metal- 
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lie phase. The intensity of the peak at the threshold, 
however, vanishes continuously at the MIT. For small 
Q < Uc/2 the opposite behavior is observed with a van- 
ishing at the transition while the integrated intensity of 
the peak remains finite. 

The behavior of a can be understood through a Fermi 
liquid analysis. The Friedel sum rule gives an exact zero- 
temperature relation between the scattering phase-shift 
and the total charge A (d-orbital plus conduction band) 
displaced when the core-hole potential is turned on^^li^^ 
namely: S = 'k/S.N /2. In the insulating phase, the valence 
band has a gap and the core-hole potential Q will not pro- 
duce a displacement of charge (AA^ = 0) unless Q is large 
enough to overcame the local repulsion (i.e. Q > U/2). 
This reflects the incompressibility of the insulator. In 
contrast, in the metallic phase, charge will gather around 
the hole to screen it even for Q <^ U . The amount of 
displaced charge reflects the vanishing compressibility of 
the metal as the MIT is approached, so that we expect 
at small Q: AA^ oc zQ. Since Fermi-liquid analysis sug- 
gests that the Nozieres and de Dominicis' value for the 
exponent [a = 2((5/7r)^] is valid throughout the metallic 
phase even in the presence of interactions, this leads to 
a parabolic behavior a oc 2 {Qz/2f - 0.05 Q'^{Uc - Uf, 
[where we used z ^ Q.?>{Uc~U)]. This expression fits our 
data remarkably well, close to Uc and for Q < Uc/2 (see 
Fig. El). 
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FIG. 4; (Color online) Core-hole photoemission spectra in the 
insulating phase {U = 3.8) for fixed Q — 2U and three values 
of ed- Left inset: in a logarithmic scale three peaks can be 
identified. Right inset: Hubbard bands in the valence spectra. 



IV. MOLECULAR ORBITAL APPROXIMATION 

The intensity and position of the peaks in the insu- 
lating phase can be understood in the molecular orbital 
approximation. The effective bath in DMFT is modeled 
with a single site having a Hubbard interaction U and 



energy e^. The resulting two site Hubbard model can 
be solved analytically. In the absence of the core-hole 
potential, the ground state is in a subspace with 2 elec- 
trons. When the core-hole potential is small {Q ^ U/2) 
the ground state is in the same charge sector and is es- 
sentially unchanged. Therefore, the XPS spectrum has a 
delta-function peak at the threshold which carries almost 
all the spectral weight (see Fig. [5]). For large Q ^ U/2 the 
final ground state is in a subspace with three electrons 
and is therefore orthogonal to the initial ground state 
((0/|0i)=0). In this regime three peaks can be identified 
in the XPS spectra and there is no peak at the thresh- 
old. The main peak has a weight 1 — — + ■ • • and 
there are two satellites at a distance ^ U + Q (shake- 
up) and ^ Q — U (shake-down). The shake-down peak 
is at a distance ed from the threshold and has a larger 
T7q 2^ + • • •) than the shake-up peak 
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FIG. 5: (Color online) Molecular orbital results for the posi- 
tion of the XPS peaks in the insulating phase. The thickness 
of the lines is proportional to the peak intensity. Parameters 
are U = 3.8D, €d = ~U/2, and = 0. 

The molecular results are in excellent agreement with 
the numerical results in the insulating phase (see Fig. [4]). 
The shake-up (-down) peak is associated to the creation 
of an electron (hole) in the upper (lower) Hubbard band 
and has an intrinsic width of order D. For a finite hole 
lifetime (finite temperature) the spectrum will have an 
additional Lorentzian (Gaussian) broadening^!. We note 
that, although it is likely to be difficult in practice, a 
measurement of the three peaks provides in principle a 
direct estimate of U. We stress that a measurement of 
the shake-down peak alone is not enough for such an 
estimation. As shown in Fig. [51 a change in ed only pro- 
duces a small redistribution of the spectral intensity and 
a global shift of the peaks relative to the threshold. The 
molecular orbital results also give a good estimation for 
the position of the peaks away from the threshold in the 
metallic phase (see Fig. [5^) . 
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FIG. 6: (Color online) Behavior of the XPS spectra at the 
edge as a function of U and Q. Left panel: coefficient a of 
the power-law divergence in the metal. Right panel: below 
the line Q ~ U /2 the spectra present a delta-function peak in 
the insulator. 

V. SUMMARY AND CONCLUSIONS 

In summary, we have studied theoretically the behav- 
ior of the core-level photoemission spectra across the 
correlation-driven MIT. Away from the photoemission 
threshold and far from the MIT, both the position and 
relative intensity of the peaks are well described by a 
molecular orbital approach. Close to the threshold or to 
the MIT, more sophisticated techniques (such as DMFT 
and NRG) are necessary to describe the spectra. The 
changes in the XPS spectra across the MIT may be used 
to detect the transition. A symmetric peak (or no peak) 
at the edge implies an insulating phase, while an asym- 
metric peak or a power-law corresponds to a metal. For 



large Q > U /2 there are three peaks in the XPS spectra. 
While the shake-up peak is usually weak and will prob- 
ably be difficult to detect, its observation in conjunction 
with the shake-down peak allows in principle for a direct 
estimation of U . 

The most interesting results concern the behavior of 
the XPS spectra at the threshold as summarized in Fig.[Bl 
In the metallic phase there is a power-law behavior with 
an exponent a — 1, where a depends on both the local 
interactions on the band U and the intensity of the core- 
hole potential Q. The power-law regime is restricted to 
an energy range ^ zD close to the threshold and disap- 
pears at the MIT. For Q < Uc/'i the exponent a depends 
strongly on U and vanishes as ~ at the transition. 
The integrated intensity, however, remains finite giving 
rise to a delta-function peak in the insulator. In con- 
trast, for Q > Uc/'i.^ a is finite at the transition but 
the integrated intensity of the peak vanishes and there 
is no delta-function peak in the insulator. Future high- 
resolution experiments might be able to test our theoret- 
ical predictions for the edge-singularity in a correlated 
metal close to the MIT. 
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